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ABSTRACT 

We report measurements of the CMB polarization power spectra from the January 2003 
Antarctic flight of Boomerang. The primary results come from six days of observation of a 
patch covering 0.22% of the sky centered near R.A. = 82.5°, Dec = —45°. The observations 
were made using four pairs of polarization sensitive bolometers operating in bands centered at 
145 GHz. Using two independent analysis pipelines, we measure a non-zero (EE) signal in the 
range 100 < £ < 1000 with a significance 4.8(t, a 2a upper limit of 8.6 fiK'^ for any (BB) 
contribution, and a 2(t upper limit of 7.0 ^K^ for the (EB) spectrum. Estimates of foreground 
intensity fluctuations and the non-detection of (BB) and (EB) signals rule out any significant 
contribution from galactic foregrounds. The results arc consistent with a ACDM cosmology 
seeded by adiabatic perturbations. We note that this is the first detection of CMB polarization 
with bolometric detectors. 

Subject headings: Cosmology, Cosmic Microwave Background, Bolometers 
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1. Introduction 

Measurements of the polarization of the Cos- 
mic Microwave Background (CMB) are a power- 
ful cosmological probe. The CMB is polarized by 
Thomson scattering (Rees 1968) during recombi- 
nation and reionization. Polarization anisotropics 
have an amplitude which is ~ 10% of the tempera- 
ture anisotropics (Bond & Efstathiou 1984). Simi- 
lar to CMB temperature anisotropics, the angular 
power spectra of CMB polarization encode cosmo- 
logical information. In recent years. CMB tem- 
perature anisotropy measurements have provided 
strong constraints on fundamental cosmological 
parameters (see e.g. Bond et al. (2003)). The 
strength of these constraints relies on the assump- 
tion that initial perturbations are adiabatic in ori- 
gin. If an admixture of isocurvature perturbations 
is allowed, then these constraints are somewhat 
weakened (Enqvist & Kurki-Suonio 2000; Buchcr 
et al. 2001). The addition of polarization informa- 
tion can constrain such isocurvature contributions 
and tighten current constraints derived from tem- 
perature anisotropics. 

With increased sensitivity, future measure- 
ments of CMB polarization will provide new inde- 
pendent constraints on the cosmological model. 
A measurement of the gravitational lensing of 
CMB polarization could provide independent con- 
straints on the neutrino mass, the dark energy 
equation of state and the nature of reionization 
(Kaphnghat et al. 2003; Hu 2002). It may also 
be possible to obtain direct evidence of inflation 
through its effect on the pattern on CMB polar- 
ization (Polnarev 1985; Crittenden et al. 1993). 

Any electromagnetic wave can be described by 
the Stokes parameters: / is the intensity, Q and 
U parameterize linear polarization, and V de- 
scribes the circular polarization. Thomson scat- 
tering does not produce circular polarization, so 
we expect V = Q for the CMB. Q and U are not 
rotationally invariant quantities. Consequently, it 
is customary to characterize CMB polarization as 
the sum of curl-free and divergence-free compo- 
nents (Zaldarriaga & Seljak 1997; Kamionkowski 
et al. 1997). Using an analogy to electromag- 
netism, the curl-free components are called E- 
modes and the divergence free components are 
called B-modes. E-modes and B-modes are re- 
lated to Q and C/ by a non-local linear transfor- 



mation. There are five observables for CMB polar- 
ization: the E-mode correlation function , {EE), 
the B-mode correlation function, {BB) , the cross- 
correlation between E-mode and B-mode polar- 
ization, (EB) and the cross-correlations between 
temperature anisotropics and polarization, (TE) 
and (TB). All of these correlations are parame- 
terized by multipole moments C/"^ where X and 
Y can represent E-modes, B-modes or tempera- 
ture anistropies. 

E-mode polarization of the CMB is primarily 
produced by scalar fluctuations on the last scat- 
tering surface, due to motion of the photon-baryon 
fluid which is induced by density fluctuations. 
However, these scalar fluctuations do not produce 
B-mode polarization on the last scattering surface. 
Tensor perturbations induced by gravity waves 
can create CMB polarization as well. Inflationary 
models generically predict a spectrum of primor- 
dial gravity waves which have an amplitude pro- 
portional to the fourth power of the energy scale 
at the time of Inflation (Turner & White 1996). 
Gravity waves produce E-mode and B-mode po- 
larization in roughly equal quantities (Seljak & 
Zaldarriaga 1997). Given current constraints on 
tensor perturbations (Seljak et al. 2005), scalar 
perturbations are expected to dominate the E- 
mode power spectrum by a factor of at least ten. 
If parity is preserved in the early universe, then 
we expect there to be no correlation between E- 
mode and B-mode polarization (i.e. {EB) = 0) 
or between temperature anisotropics and B-mode 
polarization {{TB) = 0). However, it is possi- 
ble to construct models where parity is violated 
and these correlations are non-zero (Pogosian et al. 
2002). 

In models seeded by purely adiabatic perturba- 
tions, acoustic peaks in the E-mode angular power 
spectrum should be ^ 180° out of phase with the 
acoustic peaks in the temperature anisotropy an- 
gular power spectrum, {TT). Peaks in the E-mode 
spectrum should line up with troughs in the tem- 
perature spectrum, because the scalar component 
of E-mode polarization is related to velocities and 
not densities on the last scattering surface. The 
angular spectrum of the cross-correlation between 
temperature anisotropics and E-mode polarization 
{{TE)) will show a series of acoustic peaks which 
occur between the peaks of the {TT) and {EE) 
power spectra. Measurements of the {EE) power 
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spectrum by DASI (Kovac et al. 2002; Leitch et al. 

2004) , CBI (Rcadhcad et al. 2004) and CAPMAP 
(Barkats et al. 2005) along with {TE) measure- 
ments by DASI, WMAP (Kogut et al. 2003), and 
CBI provide evidence that the assumption of adi- 
abatic perturbations is valid. 

In this paper, we report a measurement of the 
E-modc polarization power spectrum from the sec- 
ond Antarctic flight of Boomerang which took 
place in January 2003 (hereafter B03). The tele- 
scope and instrument configuration from the 1997 
test flight and the first Antarctic fiight arc dis- 
cussed in Piacentini et al. (2002) and Grill ct al. 
(2003) respectively. The instrument configuration 
for the 2003 fiight is described in Masi ct al. (2005) 
along with the data processing and the CMB 
maps. Results for the temperature anisotropy 
spectrum and the temperature-polarization cross- 
correlation are reported in Jones et al. (2005a) and 
Piacentini et al. (2005) respectively. Cosmological 
parameter constraints are reported in MacTavish 
et al. (2005). In this paper, we briefiy review the 
instrument and observations in §2. In §3, we dis- 
cuss the analysis methods used to estimate the 
polarization power spectra, and in §4 we present 
the power spectrum results. §5 provides a discus- 
sion of systematic errors, and §6 describes tests for 
foreground contamination. 

2. Instrument and Observations 

Boomerang is a balloon-borne telescope de- 
signed for long duration flights around Antarc- 
tica. In its flrst Antarctic flight (December 
1998), Boomerang measured CMB tempera- 
ture anisotropics using bolometers operating with 
bands centered at 90, 150, 220 and 410 GHz (Crill 
et al. 2003; Ruhl et al. 2003). For the 2003 flight, 
the receiver was re-designed to measure CMB tem- 
perature and polarization anisotropics with bands 
centered at 145, 245 and 345 GHz (Masi ct al. 

2005) . The results reported here come from four 
pairs of polarization sensitive bolometers (PSB's) 
operating at 145 GHz (Jones et al. 2003, 2005b). 

The 2003 flight of Boomerang was launched 
on January 6, 2003 from McMurdo Station, 
Antarctica and lasted 14 days. In this paper, 
we report on the analysis of 205 hours of CMB 
observations during the first 11 days. During this 
period 75 hours were spent scanning a large region 



(called the shallow region) comprising 3.0% of the 
sky and 125 hours on a small region (called the 
deep region) comprising 0.28% of the sky. The 
shallow region was designed to optimally measure 
the degree scale signals in (TT) and (TE), while 
deep region was designed to optimize the signal- 
to-noise ratio on the (EE) and {BB) and {EB) 
power spectra. For these spectra almost all the 
statistical weight comes from the deep region. 

Figure 1 shows the distribution of integration 
time over the observation region and the sky cuts 
used in the analysis. The time per pixel in the deep 
region is a factor of ^ 20 longer than in the shallow 
region. For the spectral analysis, we use 1.8% of 
the sky for the shallow region and 0.22% of the 
sky for the deep region. These choices were made 
so that the coverage was roughly uniform in time 
per sky pixel and for the different channels in the 
focal plane. Another consideration was that the 
deep and shallow observations could each be split 
in half (in time) and still cover their respective sky 
cuts (which is useful for systematic tests). 

3. Data Analysis Methods 

We used two independent pipelines for the map- 
making and polarization power spectrum estima- 
tion. Masi et al. (2005) describes the bulk of the 
data analysis from raw data to CMB maps includ- 
ing: raw data cleaning, detector characterization, 
pointing reconstruction, calibration, beam mea- 
surement, noise estimation and making polarized 
maps. In this paper, we limit the discussion to as- 
pects of power spectrum estimation for polarized 
datasets. 

The two pipelines are independent to a high 
degree. In the following, "NA pipeline" refers to 
the pipeline primarily based in North America and 
"IT pipeline" refers to the pipeline developed in 
Italy. Each team had many choices to make about 
instrument characteristics and data analysis tech- 
niques. The boundaries of the shallow and deep 
region sky cuts used by both teams are identical, 
but differences in the data cleaning causes slight 
differences in integration time. As we will show in 
the §4, the two pipelines yield compatible answers; 
this is a testament to the robustness of the data 
set. 
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Fig. 1. — Sky coverage from the 2003 Boomerang flight. The top panel shows the large region covered during the 
first part of the flight (the shallow region) and the bottom panel is the smaller region covered during the second half 
of the flight (the deep region). In the top panel, the outer set of black lines shows the sky cut used in the shallow 
mask. The inner outline shows the outline of the deep region sky cut. The shallow scans covered 3.0% of the sky, 
and sky cut used for the CMB analysis covers 1.8%. The deep region observations covered 0.28% of the sky and the 
outlined region covering 0.22% of the sky was used for the CMB analysis. The integration time per pixel for the deep 
observations is roughly 20 times longer than integration time per pixel during shallow observations. In both panels, 
the small circles represent regions of map which are excised due to the presence of known point sources. 
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3.1. Power Spectrum Estimation 

Both pipelines are polarized extensions of the 
Monte Carlo based MASTER method (Hivon 
et al. 2002) first used on B98 (Netterfield et al. 
2002). These techniques rely on spherical har- 
monic transformations done on a partial map of 
the sky. For polarization data, the Q and U maps 
are expanded as a function of spin- 2 spherical 
harmonics 

(g ± iU){h) = ^(af„ ± tafj±2Yim{n), (1) 

Im 

where and are the coefficients for E-mode 
and B-mode polarization respectively. These co- 
efficients can be calculated in a manner similar to 
Legendre transformations, 



dnW{h) {Q + tU){h)+2Yi„,{n) 



+ {Q-iU){n)^2Yim{n) 
1 



(2) 



{Q-iU)ih)^2Yim{n) 
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where VF(n) is an arbitrary weighting function and 
the integral extends only over the observed portion 
of the sky. From these transforms, we can build 
three observablcs: 
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where the E-mode power spectrum, Cf^ 

the B-modc power spectrum, and C™ the cross- 
correlation between E-mode and B-mode polariza- 
tion. C™ is expected to be zero if parity is pre- 
served in the early universe. Our estimates of the 
cross-correlations between temperature and polar- 
ization [Cj'^ and Cj'^) arc discussed in Piacentini 
et al. (2005). 

For spherical harmonic transforms done on the 
cut sky the measure of Ci is biased; we describe 
them as pseudo-C^'s {Ci). For the polarization 



power spectra, the relationships between full-sky 
Ci and Cf are expressed as 
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where represents the full-sky power spec- 

trum. Eg is the beam window function, is the 
transfer function measured by signal-only Monte 
Carlo simulations, -/V^^ is the noise bias measured 
by noise-only Monte Carlo simulations, +Kiii is 
the primary coupling kernel and -Kw describes 
the geometric leakage between E-modes and B- 
modes (Chon et al. 2004). Both pipelines use 
roughly 500 Monte Carlo simulations of signal- 
only and noise-only data streams to estimate the 
signal transfer function and noise bias respectively. 
A similar number of signal-|-noise simulations can 
be used to estimate the uncertainty on the spectral 
estimate and check for bias in the pipeline. 

Since we observe a small portion of the sky, we 
are not able to measure individual multipole mo- 
ments. Instead, we parameterize the power spec- 
trum as a piecewise continuous function 



— lb ! 



(10) 



where q^"^ is the bandpower deviation over a 
range (A^)f, and c\^^^^ is a shape parameter. 
Common choices for the shape parameter are those 
that keep Cj'^'' constant over the band, those that 
keep -f- 1)C^^^^ / {27:) constant over the band 
(i.e. the flattened spectrum) or those that rep- 
resent a theoretically motivated power spectrum 
(e.g. ACDM concordance model). The choice of 
parameterization depends in part on the nature of 
the expected signal and the noise in the maps. 
The output bandpower C^'^ [Ci = £(£ + 
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l)Cf/27r) is then a function of and 
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When comparing to a model, the expected band- 
power deviation can be written as 



(12) 



where Wj^ is the bandpower window function, 
' = e{e+l)cf^ /2'K andl[fi] is the logarithmic 
integral (Bond et al. 2000) 
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For general shape functions, we get 

{qb) = 



(13) 
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and we can recover (Cb) using equation 11. If 



C 



constant, then we have 



{Cb 



3.2. NA Pipeline 
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For the NA pipeline, a quadratic estimator 
(Bond et al. 1998) is used to iteratively solve for 
the bandpowcrs and their uncertainty. This esti- 
mator (called Xfaster) is capable of solving for the 
power spectra of a single map or any combination 
of two or more maps (which can be overlapping) 
while accounting for all correlations between those 
maps (Contaldi et al. 2005). 

For BOS, Xfaster is used to solve for the com- 
bined power spectrum of the shallow and deep re- 
gion data (the combined power spectra are called 
the 2Mask spectra) . Separate maps are made from 
the shallow and deep region observations; this in- 
sures that the only correlations between the maps 
are due to sky signal. When we perform the spher- 
ical harmonic transformations, we use a uniform 
pixel weighting for the shallow map, and pixels in 
the deep region are weighted by the inverse square 

— 1/2 

root of noise in that pixel {(Jp^J )■ An effective 



noise weighting is also applied in the spectrum es- 
timation process due to noise bias of each map 
(Nf-^). This is an efficient way to account for the 
imbalance of integration time per pixel between 
the shallow and deep maps. 

Although the shallow region docs not con- 
tribute much statistical weight to the polariza- 
tion spectra, it does significantly reduce the sam- 
ple variance in the (TT) and (TE) spectra. The 
Xfaster method is used for the polarization spec- 
tra so that we can derive consistent correlation 
matrices between all spectra for use in parameter 
estimation (MacTavish et al. 2005). 

3.3. IT Pipeline 

To solve for the power spectrum, the IT pipeline 
uses a method similar to that described in Hivon 
et al. (2002) but adapted for polarization spec- 
tra (see e.g. Kogut et al. (2003) and Challinor 
& Chon (2005)). Errors bars and correlation ma- 
trices are calculated using bandpowers resulting 
from signal-|-noise Monte Carlo simulations. The 
IT pipeline performs the (EE), (BB) and (EB) 
analysis on the deep region maps with the pixels 
weighted by the inverse square root of the noise 
i'^pix )• The IT maps also include data from shal- 
low observations which fall inside the deep region. 

3.4. Testing Goodness-of-Fit 

To test how well particular models fit our data, 
we use the standard likelihood ratio technique in 
a manner similar to recent work done by DASI 
(Kovac et al. 2002; Leitch et al. 2004) and CBI 
(Readhead et al. 2004). Specifically, we calculate 
the logarithm of the ratio of the peak of the likeli- 
hood to the likelihood of a model Ai parametrized 
by bandpowcrs Cf,: 



A(C6 



In 



LjCb) 
L{Cb) 



(16) 



where Cb are the maximum likelihood bandpow- 
ers. The larger this ratio is, the worse the model 
A4 fits the data. In this paper, we are primarily 
concerned with comparisons to the null hypoth- 
esis {Cb = for all b) and our fiducial ACDM 
model (the best fit to the WMAP (TT) spectra 
from Spergel et al. (2003)). 

In the approximation that the likelihood func- 
tion L{Cb) is a multivariate Gaussian near its peak 
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(which is generally a good approximation when the 
signal to noise ratio in a band is ^ 1), we have 
A = If we assume that a given model M 

is true, the probability of observing A exceeding a 
particular value of is given by the "probability to 
exceed" (PTE) 

where N is the number of parameters and T(x) 
is the complete gamma function. For example, if 
PTE=5%, then we can reject the hypothesis that 
the model A4 is true with 95% confidence. 

Since our power spectrum estimators do not cal- 
culate full likelihood values, we use the offset log- 
normal function, = ln{Cb + Xb), to approximate 
the likelihood function (Bond ct al. 2000) where Cb 
is the bandpowcr and Xb is the offset parameter. 
The likelihood is calculated by 

(7b=ACb/{Cb + Xb) (18) 
AZb = ln(Cfc + Xb) - ln(Cb + Xb) (19) 
-21ni(Cfc) ^Y.^^h<yb'Gbb'(Tb''^Z^' (20) 

bb' 

where Cb is the maximum likelihood bandpower 
and Gbb' is the bandpowcr correlation matrix 
which is normalized to unity on the diagonals. In 
this parameterization, the likelihood is normalized 
to the peak value (i.e. lnL(C;,) = 0). 

In the results to follow, we compute A and PTE 
separately for the (EE) , (BB) and (EB) spectra. 
In other words, when performing this test on one 
spectrum, we marginalize over the other spectra 
by excising the correlations between spectra from 
the inverse Fisher matrix. 

4. Results 

4.1. Narrow Band Analysis 

The (EE), (BB) and (EB) power spectra re- 
sults from both pipelines are shown in Figure 2 
and listed in Table 1. The multipole range shown 
is from 101 < £ < 1100; information on scales 2 < 
£ < 100 and £ > 1101 is discarded. These data, 
XbS, window functions and correlation matrices 
are available at http://cmb.case.edu/boomerang 
and http://oberon.romal.infn.it/boomerang. 

Although we use a shape function which is con- 
stant in £{£ + l)Cf/(27r), wc choose to plot the 



results in terms of {2£ + l)C£/(47r) to emphasize 
our sensitivity in the range 100 < £ < 1100 and 
relate the power spectrum directly to the r.m.s. 
CMB signal 

^cA/B = X! 4^ '^i- (21) 

I 

The results from the two pipelines show a high 
degree of consistency. Additionally, the NA deep- 
only power spectra are nearly identical to the 
2Mask spectra plotted here. Different choices in 
the data processing (e.g. time domain filtering, 
the number of map iterations, detector relative 
calibrations) could lead to the small discrepancy 
between the NA and IT results. 

In Figure 2, only the {EE) spectrum appears 
to be significantly different from zero. To quan- 
tify this, we calculate the A statistic for the as- 
sumption of zero polarized signal. The results in 
Table 2 show that the (EE) result is inconsistent 
with zero signal and that {BB) and (EB) are both 
consistent with zero. Similarly, we compare (EE) 
and (BB) to the E-mode power spectrum from 
the fiducial ACDM model. In this case, the (EE) 
result is in good agreement with the fiducial E- 
mode spectrum while (BB) is not. For these cal- 
culations, the IT data and the NA (EB) data are 
taken to have Xf, = 0. 

4.2. Wide Band Analysis 

To assess the raw significance of our (EE) re- 
sult and to set an upper limit for {BB) and {EB) 
spectra, we perform a wide band analysis over a 
range 201 < £ < 1000 (the actual analysis uses 
three bins with bins defined by 2 < € < 200 and 
1001 < £ < 1999 used as "junk" bins). We used 
four different shape functions for the bandpowers: 
constant in C(, constant in {2£+l)Ce/4:Tr, constant 
in £{£ + l)C(-/27r and the fiducial ACDM model. 
Sky cuts and spectrum estimation are the same as 
those used in the narrow band analysis. 

Table 3 shows the bandpowcr results for all 
cases of the wide band analysis and Table 4 shows 
the A statistic and PTE for each case. In all cases, 
the {BB) and {EB) results are consistent with 
zero signal while the {EE) signal is significantly 
non-zero. The choice of shape function and the 
fine details of the bandpower estimator have an 
effect on the output bandpowers. The NA and 
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Fig. 2. — Angular Power Spectra Results. From top to bottom the panels show the (EE), (BB) and (EB) power 
spectrum results from the NA (blue circles) and IT (red squares) pipelines. The solid line in the (EE) plot is the 
best fit ACDM model to the WMAP {TT) resuhs (Spergel et al. 2003) and the dashed hue in aU plots represent 
zero-signal. From these plots and the statistical tests in Table 2, it is clear the (BB) and (EB) are consistent with 
zero signal while {EE) is consistent with ACDM, but inconsistent with zero signal. The ACDM (EE) model predicted 
by the WMAP (TT) results is nearly identical the best fit {EE) model predicted from recent (TT) results including 
BOS (MacTavish et al. 2005). 



IT bandpowcrs agree closely with the results from 
Monte Carlo simulations of each method using the 
fiducial ACDM model as the input. When we pa- 
rameterize the spectra as flat in £{£ + l)C^/(27r), 
the simulations also show that the difference in er- 
ror bars is consistent with the differences between 
the NA and IT estimators. This difference is illus- 
trated by the top panel of Figure 3 which shows 
the wide-band window function used by the NA 
Xfaster estimator. Xfaster performs an effective 
Wiener filter where most of the statistical weight 



comes from the lower edge of the band. The IT 
estimator uses a fiat window function which leads 
to a less optimal result. The bottom panel of Fig- 
ure 3 shows how the NA window functions varies 
with the choice of shape function. 

For each parameterization, the A statistic and 
PTE results from the wide band analysis are re- 
ported in Table 4. For a single bandpower, the 
significance of the detection can be calculated by 
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-1.0 


3.8 


{EB) 


301 


400 


350 


-7.6 


4.8 


-6.8 


5.0 


401 


500 


450 


-6.6 


6.3 


-5.4 


6.6 


501 


700 


600 


-1.1 


7.6 


-6.9 


6.9 




701 


900 


800 


-9.1 


19.1 


17.9 


16.2 




901 


1000 


1000 


122.1 


50.4 


8.2 


44.3 



Table 1: Power spectra results for {EE}, {BB) and {EB) from the NA and IT analysis pipelines. The first three 
columns define the i'-bins used and the next four columns state the bandpowers and errors from both pipelines. Cj, 
and AC(, {Ci = + 1)Ci/2-r) have units ^K^ . Power spectra, a;i,'s, window functions and correlation matrices are 
available at http://cmb.case.edu/boomerang and http://oberon.romal.infn.it/boomerang. 



where S is the detection significance quoted in 
units of (7. With the NA results, we find that the 
shape function which is constant in Cil2-K pro- 
duces the highest significance, but its significance 
is only slightly higher than what we find when us- 
ing {2i+l)Ci/A'n or ACDM as the shape function. 
For the case parametrized as constant in Cf, we 
find that the {EE) bandpower is consistent with 
a 4.8(7 detection. For this same case, we quote 2ct 
upper hmits of 8.6 ^iK'^ for {BB) and 7.0 iiK'^ for 
{EB). 

5. Systematic Errors 

Given the small amplitude of the polarization 
signal, we need tight control on systematic errors. 
For most systematic errors, we would expect them 



to contribute equally to {EE) and {BB) which 
could also lead to a non-zero {EB). The fact 
that {BB) and {EB) are consistent with zero gives 
credibility to the {EE) result. To further estab- 
lish the robustness of our result, we performed 
two types of internal consistency checks and a 
suite of Monte Carlo simulations to determine lim- 
its on systematic errors due to instrument mis- 
char actcrization. 

5.1. Internal Consistency Tests 

To check the consistency of our result, we per- 
formed jackknife tests which are done by splitting 
the data in half, making maps Ai and A2 from 
each half and measuring the power spectrum of 
(Ai — A2)/2. If this power spectrum is consis- 
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Fig. 3. — Window functions from the NA wide band results in Tables 3 and 4. In the top panel the solid blue 
line is the {EE) window function for the band 201 > I < 1000. The dashed red line characterizes the leakage of 
E-modes into B-modes [E B) which has a maximum value of ~ 0.01. The {BB) window function and B ^ E 
leakage window function are similar to those plotted here. The low amplitude of the E —> B and B —> E shows 
that BOS is able to separate E-mode and B-mode polarization. In the bottom panel, the window functions are shown 
for the different power spectrum parameterizations (i.e. C^^') used in the bandpower estimation. The shape of the 
(EE) window function indicates the effective weight applied to each multipole moment. For all cases used in the NA 
analysis, the window function is significantly different than the flat band used in the IT wide band analysis. This is 
due to an effective Wiener filter which weights each multipole by Cf^jiCi + A'^)'^ where Ni is the noise at a given 
multipole, C\ is the shape function and Ci is the expected signal. In a given band the expected signal depends on 
the form of C^"^'. 



tent with zero then the dataset is considered to 
be internaUy consistent. We performed two sets 
of jackknife tests. The first test involves splitting 
the data in time (caUed the (h\ — h2)/2 test). The 
second test is done by comparing detectors on the 
left and right side of the focal plane (called the 



(WX-YZ)/2 test). The [hi - h2)/2 is sensitive to 
time- varying systematic problems while the (WX- 
YZ) /2 is sensitive to problems affecting individual 
channels. 

The mapmaking process for the {hl — h2)/2 test 
is different in the NA and IT pipelines. With the 
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NA pipeline, we make maps from the first and 
second fialf of the shallow observations (A^j^ and 
A^2)j Eind first and second half of the deep obser- 
vations (A^j and A^2)- We then use Xfaster to 
estimate the combined power spectrum of (A^-^ — 
Af2)/2 and (A^^ - Af2)/2. With the IT pipeline, 
we make a combined map from the first half of the 
shallow observations and the first half of deep ob- 
servations A^f . Similarly, A^^ is made from sec- 
ond halves of the shallow and deep observations. 
The IT esimator is then used to estimate the power 
spectrum of (A^f — Af^)/2 on the deep region 
mask. 

As discussed in Masi et al. (2005), each side 
of the focal plane has two PSB pairs at 145 GHz 
which were oriented so that the left and right sides 
of the focal plane could measure Stokes Q and U 
independently. The (WX-YZ)/2 is done by taking 
the difference of maps made from from the left 
(WX) and right (YZ) sides of the focal plane. In 
the same manner used for the {hi — h2)/2 test, 
the NA pipeline makes separate (WX-YZ)/2 maps 
from the shallow and deep observations while the 
IT pipeline makes a combined (WX-YZ)/2 map 
from the deep and shallow data. For each pipeline, 
the power spectra are estimated in the same way 
as in the {hi — /i2)/2 test. 

Figure 4 shows the results for the {hi — h2)/2 
and the (WX-YZ)/2 tests, and Table 5 shows 
and PTE's calculated from those results . For each 
pipeline, both jackknife tests are consistent with 
zero for all three spectra. These tests put strong 
limits on systematic problems. 

5.2. Simulation of Instrument Characteri- 
zation Errors 

Mis-characterization of instrumental parame- 
ters is a potential source of systematic error in the 
power spectra. The primary parameters of con- 
cern are: beam size, calibration, polarization effi- 
ciency, detector time constant and polarization an- 
gle. An error in beam size leads to a bin-dependent 
scaling factor. Errors in the absolute calibration 
and polarization efficiency lead to an overall scal- 
ing factor. Errors in relative calibrations or detec- 
tor time constants lead to leakage of CMB tem- 
perature anisotropics into the polarization signal. 
An error in the polarization angle mixes the Q and 
U Stokes parameters. 



The measurement of B03 instrument parame- 
ters is described in Masi et al. (2005) and the un- 
certainties on those parameters are shown here in 
Table 6. To estimate the error induced by poten- 
tial errors in relative calibration, time constant, 
polarization efficiency and polarization angle, we 
performed a suite of signal-only Monte Carlo sim- 
ulations. For each parameter, we performed 145 
simulations starting with the same simulated sky 
map (a realization of the fiducial ACDM model). 
We then create a time-ordered datastrcam where 
the value of the parameter is randomly varied for 
each detector. The values are drawn from a dis- 
tribution representing our uncertainty on that pa- 
rameter. We then analyze this data stream using 
the measured parameter values. For each Monte 
Carlo, we estimate the power spectrum using a 
technique similar to that used by the IT pipeline. 
We then compute the systematic error bar by tak- 
ing the standard deviation of Monte Carlo results. 
Figure 5 shows the results of the simulations. The 
induced systematic error bars are less than 10% of 
the bandpower uncertainty. On most scales, the 
polarization angle is the dominant source of error. 

Although uncertainty in the beam size is a rel- 
atively benign problem, beam differences between 
elements in a PSB pair and structure in the cross- 
polar beam pattern of a given detector could lead 
to irreducible leakage of temperature anisotropics 
into the polarization maps. Given the off- axis 
structure of the Boomerang optics (where the 
axis of symmetry is vertical) , the beam mismatch 
between elements in a PSB pair depends on the 
polarization angles of each element. For example, 
detectors oriented at 45° and —45° with respect 
to horizon have well matched beams while detec- 
tors oriented at 0° and 90° will have slightly differ- 
ent beams. A physical optics simulation confirmed 
that the latter case has the worst mismatch among 
all four PSB pairs. For this worst case, we calcu- 
late the differential beam window function and es- 
timate the leakage of temperature anisotropy into 
polarization using the fiducial ACDM model. We 
find this signal to be smaller than our measured 
{EE) signal by a factor of 10^ (or more) on the an- 
gular scales that we are sensitive to. This calcula- 
tion represents the worst case scenario because sky 
rotation should reduce this contamination some- 
what. 

Differences between the cross-polar beams 
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Fig. 4. — Results of Jackknife Tests. The left side shows the results for the {hi — h2)/2 test and the right side shows 
the results for the (WX-YZ)/2 test. The blue circles are results from the NA pipeline and the red squares are results 
from the IT pipeline. Table 5 shows the and PTE calculated from these results. For both tests, all three spectra 
are consistent with zero signal. 



in a FSB pair could also lead to temperature 
anisotropy leakage. The integrated cross-polar 
beam for a given detector is a factor of ~ 200 
smaller than the integrated co-polar beam. Here 
if we take the worst case scenario (i.e. that the 
throughput of cross-polar beam difference is twice 
the throughput of the cross-polar beam of one de- 
tector), a naive estimate of the (TT) leakage into 
{EE) gives {EE)i^^f^ < 10-*{TT) which is < 1% 
of the observed (EE). 



6. Foregrounds 

Folarized emission from galactic and extra- 
galactic sources are another potential source of 
contamination (dc Olivcira-Costa 2004; Tucci 
et al. 2005). Currently, not much is known 
about diffuse polarized emission in the frequency 
range 30 < ^ < 400 GHz. Synchrotron emis- 
sion is expected to be highly polarized, but the 
power steeply decreases with increasing frequency 
(a ~ —3). Recent observations with the ACTA 
telescope have detected polarized synchrotron 
emission in a small patch near the edge of our 
deep region at frequencies of 1.4 GHz (Bernard! 
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Fig. 5. — Propagation of measurement errors in in- 
strumental parameters to (EE) error. The hatched 
bands show the upper edges of the (EE) error (due 
to noise) for the bands in the multi-bin NA results 
(Figure 2). The other symbols show errors in rela- 
tive calibration (stars), bolometer time constant (up- 
ward triangles) and polarization angle (downward tri- 
angles). Because errors in absolute calibration and 
polarization efficiency are multiplicative factors which 
act identically on each bin, their effect is left off this 
plot and reported instead in Table 6. 

et al. 2003) and 2.3 GHz (Carretti et al. 2005). 
A naive extrapolation of these synchrotron results 
to 145 GHz predicts a signal of 0.2 r.m.s. 
compared to a ~ 3 ^lK r.m.s. expected from the 
fiducial ACDM model given our beam. 

From starlight polarization measurements (Fos- 
alba et al. 2002), dust is expected to be less than 
10% polarized, with a spectral index a ^ 1.7, 
but it could be higher depending on the nature of 
the galactic magnetic field (Wright 1987). Results 
from Archeops (Benoit et al. 2004; Ponthieu et al. 
2005) measure a polarization fraction of 5 — 10% 
for dust clouds near the galactic plane, but the 
results are not sensitive enough to place strong 
limits on degree scale dust polarization away from 
the galactic plane. 

From the 145 GHz data alone, we are highly 
confident that our E-mode polarization signal is 
dominated by the CMB and not foreground emis- 
sion. Foreground emission should produce nearly 



equal parts E-mode and B-mode polarization. 
Non-detection of any B-mode signal (Figure 2 
and Tables 1-4) strongly implies a lack of fore- 
ground polarization. In Piacentini et al. (2005), 
the B03 {TB) signal is consistent with zero while 
{TE) is consistent with ACDM. Further evidence 
is obtained by cross-correlating an IRAS dust in- 
tensity map with the 145 GHz polarization data. 
Both the (TiRAsBBOi) and {TirasEbos) are con- 
sistent with zero. We find that the consistency of 
the deep-only and combined shallow-|-deep power 
spectra rule out large dust polarization signals in 
regions nearer to the galactic plane. 

In Masi et al. (2005), we characterize the dust 
emission by comparing the three B03 intensity 
maps to dust templates from Schlegel et al. (1998). 
In the deep region, we detect a dust intensity cor- 
relation at 2.5ct with our 345 GHz channels, but 
find an upper limit of 4 I-iKcmb r.m.s. for dust in- 
tensity at 145 GHz. If we assume that dust is 10% 
polarized, we get an upper limit of 0.4 ^Kqmb 
r.m.s. for the dust polarization signal at 145 GHz. 
A polarized analysis of 245 and 345 GHz data will 
be discussed in a future work. 

Lastly, the difference of the intensity maps at 
145 and 345 GHz shows what appears to be three 
small regions of diffuse dust emission which appear 
to be correlated with IRAS emission (see Figure 
27 of Masi et al. (2005)). As one final test, we 
perform a spectrum analysis on a sky cut where 
we excised square blocks centered on these clouds. 
The boundaries of the blocks are reported in Ta- 
ble 7. The resulting polarized power spectra are 
identical to the ones reported here. 

7. Conclusions 

In this paper, we report the measurement of the 
{EE) {BE) and {EB) polarization power spectra 
from the 2003 flight of Boomerang. Consistent 
results have been obtained from two different data 
analysis pipelines. These results have passed a 
wide variety of systematic tests and the induced 
error from instrumental uncertainties is negligi- 
ble. The {BB) and {EE) results are consistent 
with zero signal, as expected in ACDM models 
dominated by scalar adiabatic perturbations. The 
{EE) results are consistent with existing measure- 
ments (Figure 6) and a good fit to the {EE) signal 
expected from the ACDM model which is the best 
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fit to the WMAP (TT) results. Several tests using 
higher frequency channels and dust maps, in addi- 
tion to the fact that (BB) and (EB) arc consistent 
with zero, argue that it is very unlikely that the 
(EE) result is contaminated by galactic emission. 
This detection of (EE) is the first by a bolometric 
polarimeter, and thus bodes well for the future of 
CMB polarimetry using bolometric detectors. 
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Fiducial Model 


No Polarization 


NA 


IT 


NA 


IT 


A PTE 
{EE) 2.6 0.63 
{BB) 6.6 0.07 
{EB) - 


A PTE 
3.2 0.49 
9.2 0.01 


A PTE 
11.8 1.3 X 10-3 
1.3 0.92 
5.8 0.11 


A PTE 
7.2 0.05 
2.0 0.78 
2.6 0.64 



Table 2: Table of values for the A statistic and "probability to exceed" (PTE) calculated from the narrow band 
results in Figure 2 and Table 1. We compare {EE) and {BB) to the fiducial ACDM model (the best fit to the 
WMAP {TT) from Spergel et al. (2003)) and we test the null hypothesis (no polarization signal) on all three spectra. 
Both data sets are assumed to have 7 degrees of freedom. For both data sets {EE) is a good fit to fiducial model 
and inconsistent with the null hypothesis. {BB) is not a good fit to fiducial model, but is consistent with no signal. 
{EB) is consistent with zero signal. 



Shape Pipeline (Cf^) 


{EE) 
Cb ACb 


{BB) 
Cb ACb 


{EB) 
Cb ACb 


1{1+1)Ci/{2tt) NA 2mask 9.94 
IT deep 15.5 


11.5 3.0 
14.2 6.3 


2.0 2.2 
8.0 5.7 


-4.00 1.9 
-0.2 3.9 


C£ NA 2niask 19.4 


23.4 5.2 


3.3 4.3 


-4.7 3.5 


(2^+l)C£/47r NA 2niask 14.9 


17.5 4.0 


2.5 3.15 


-5.1 2.6 


ACDM NA 2niask 15.4 


16.4 3.8 


3.3 3.0 





Table 3: Wide band analysis results for 200 < £ < 1000. The Shape column refers to the band power parameterization: 
flat in £{£+l)Ce/{2Tr), flat in Ce and ACDM. The column labeled (Cjf ^) is the {EE) expectation value for an ensemble 
of Monte Carlo simulations using the BOS instrument noise and the fiducial ACDM model as the CMB input. {C^^) , 
Cb and ACf, have units fj.K'^ and are quoted in terms of £{£ + l)C</(27r). {EB) is not calculated for the ACDM shape, 
since it can go positive and negative. 



Shape Mask 


{EE) 
A PTE 


{BB) 
A PTE 


{EB) 
A PTE 


i{e + l)Q/(27r) NA 2mask 
IT deep 


8.6 3.4 X 10-^ 
2.5 0.02 


0.42 0.36 
1.0 0.16 


2.2 0.04 
1.0 X 10-3 0.97 


Ci; NA 2mask 


11.4 1.9x10-^ 


0.30 0.44 


0.94 0.17 


(2£+l)C£/47r NA 2mask 


10.7 3.7x10-'' 


0.32 0.42 


1.9 0.05 


ACDM NA 2mask 


10.4 5.3 X 10-*^ 


0.61 0.27 





Table 4: Values of the A statistic and the PTE from the comparison of the wide band the bandpowers in Table 3 to 
a model with zero polarization signal. The {EE) data is inconsistent with zero, while {BB) and {EB) are consistent 
with no polarization signal. 



16 







NA 




IT 


Spectrum 




PTE 




PTE 


{EE) {hi - h2)/2 


6.9 


0.44 


10.6 


0.16 


{BB) {hi - h2)/2 


3.0 


0.89 


6.0 


0.54 


{EB) {hi - /i2)/2 


3.0 


0.89 


6.9 


0.44 


{EE) (WX-YZ)/2 


5.5 


0.60 


5.9 


0.75 


{BB) (WX-YZ)/2 


6.7 


0.46 


1.5 


0.997 


{EB) (WX-YZ)/2 


5.1 


0.65 


6.0 


0.74 



Table 5: Reduced and "probabilty to exceed" calculated from the jackknife results (Figure 4). Both 
pipelines have 7 degrees of freedom. 



Parameter Uncertainty Induced Error on Spectrum 

beam FWHM (TW 2.5% at £ = 500, 10% at £ = 1000 

absolute calibration 1.8% 3.6% 

polarization efficiency 3% 4% 

relative calibration 0.8% See Figure 5 

polarization angle 2° See Figure 5 

time constant 10% See Figure 5 



Table 6: Instrument parameters, the uncertainty on their characterization and the induced error on the {EE) power 
spectrum. Errors in beam size, absolute calibration and polarization efficiency result in a re-scaling of the power 
spectrum. Errors in relative calibration, polarization angle and time constant are more complicated; see Figure 5 for 
the induced error bars. 



R.A. 


limits 


Dec. 


limits 


84° 


85° 


-48.5° 


-47.25° 


87.5° 


88.75° 


-48.5° 


-47.25° 


87.5° 


88.5° 


-49.5° 


-50.5° 



Table 7: Regions of potential dust contamination where found by taking the difference between the B03 145 and 
345 GHz intensity maps. The polarization spectra were re-calculated with the data in these regions excised. The 
resulting spectra are identical to those in Figure 2 and Table 1. 
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